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We investigate multiexciton bound states in a semiconducting phase of divalent hexaborides. 
Due to three degenerate valleys in both the conduction and valence bands the binding energy of a 
6-exciton molecule is greatly enhanced by the shell effect. The ground state energies of multiexciton 
molecules are calculated using the density functional formalism. We also show that charged impuri- 
ties stabilize multiexciton complexes leading to condensation of localized excitons. These complexes 
can act as nucleation centers of local moments. 
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I. INTRODUCTION 

The recent discovery of unusual high-temperature 
weak ferromagnetism in lightly doped divalent hexa- 
borides, CaBg, SrBg, BaB@, 0] has raised the need to un- 
derstand the properties of these compounds in a broader 
perspective. The unique features of the novel type of 
ferromagnetism in Cai-jXa^Bg are (i) there are no par- 
tially hlled d- or /-orbitals, (ii) magnetism appears only 
at finite doping < x < 2%, and (iii) tiny magnetic mo- 
ments (~ per doped electron) are very robust and 
develop at temperatures as high as 600-1000 K. All these 
features are consistently explained by ferromagnetism of 
a doped excitonic insulator O). Thus, not only is the 
magnetic phase quite unusual, but also the undoped sto- 
ichiometric hexaborides may also exhibit a novel type 
of ground state — a condensate of bound electron-hole 
pairs or excitons. The theory of excitonic insulators has 
been developed long time ago, see ||. However, divalent 
hexaborides seem to be the first experimental realization 
of an excitonic insulator made possible by their unique 
band structure. Band structure calculations [|||| predict 
a small overlap of valence and conduction bands at the 
3 symmetry related nonequivalent X points in the cubic 
Brillouin zone. There is a certain ambiguity both from 
theoretical and experimental ||:@] points of view re- 
garding whether the stoichiometric divalent hexaborides 
have a small band overlap or a small band gap. Precise 
treatment of correlation effects in electron-hole plasma 
beyond the standard approximations of band struc- 
ture calculations predicts a first order transition from a 
semiconductor to a semimetal under pressure. Therefore, 
stoichiometric hexaborides may be on a semiconducting 
side of the phase diagram close to the first order metal- 
insulator transition. 

Motivated by the unique band structure with equal 
number of degenerate valleys in conduction and valence 
bands with different symmetries, we investigate here the 
formation of multiexciton molecules and complexes in 
these materials, assuming that there is a finite band gap 
between conduction and valence bands. In the effective 



mass approximation the valley index appears as an ex- 
tra quantum number. By analogy with nuclei we ex- 
pect a shell structure of the single particle orbits where 
six electrons and six holes with different spin and valley 
quantum numbers can occupy the lowest shell. Hence, 
the 6-exciton molecule is expected to be the most stable 
configuration. 

Binding energy calculations of several excitons have 
been reported in the late seventies for multiexciton com- 
plexes bound to impurities in Si and Ge § ||. It was 
recognized that Hartree-Fock calculations do not work for 
multiexciton complexes because the main contribution to 
binding energy comes from correlation effects and that 
the density functional formalism provides a useful tool 
to include such correlation effects. In this paper we con- 
sider both multiexciton molecules and multiexciton com- 
plexes bound to impurities and calculate their binding 
energies within the framework of the density functional 
theory. The Kohn-Sham equations are single particle 
Schrodinger equations for electrons and holes moving in 
a self-consistent potential. These equations connect the 
density functional formalism to the shell model of multi- 
exciton complexes and molecules ||. 

The binding energies of exciton molecules provide 
information on whether a Bose gas of multiexciton 
molecules is stabilized between semiconducting and 
semimetallic phases when the band gap is varied from 
positive to negative values. In the previous study by 
two of us H an intermediate phase of a free exciton gas 
stabilized by intervalley scattering processes was consid- 
ered. Binding of several excitons is an alternative mech- 
anism for appearance of an intermediate phase between 
a semiconductor and a semimetal. The binding energy 
per electron-hole pair in the molecule has to be compared 
to the minimal ground state energy of the electron-hole 
liquid, which according to Ref. M is S™" 1 = —1.5525a, 
in units of the excitonic Rydberg E x . As another ap- 
plication of our results we suggest that the formation of 
multiexciton complexes near donor impurities (e.g., La- 
substitutions) can be a source of local magnetic moments, 
which appear in a semiconducting state. Such moments 
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could explain the unusual NMR relaxation rate measure- 
ments in the hexaborides [fL3|. 



II. METHOD 

We use the effective mass approximation consider- 
ing electrons near minima in the conduction band and 
holes near the maxima of the valence band as oppo- 
sitely charged quasiparticles with anisotropic effective 
masses which interact via a screened Coulomb poten- 
tial V(r) — e^/er, e being the static dielectric con- 
stant. For CaBg the principal values of the effective 
mass tensor in units of the bare electron mass are given 
by ml = 0.504, = 0.212 (conduction band) and 
m\ — 2.17, — 0.206 (valence band). The model 
Hamiltonian is 

H = f + U , 
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where A = (k, v, a) with k = e, h; v = 1, 2, 3 is the valley 
index; a = ± is the spin index; and q\ is +e for positively 
charged holes and — e for negatively charged electrons. 
The kinetic energy operator K\ for one valley is 
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The expressions for the valleys 2 and 3 are obtained 
by cyclic permutations of x, y and z. We investigate 
the ground state energy of this Hamiltonian for a given 
number of electrons and holes. To treat this problem 
we use the density functional formalism which is pre- 
sented in the next section. In the following all equations 
are written in dimensionless variables. A natural set of 
units is the excitonic Rydberg for the energy scale E x = 
/ie 4 /2h 2 e 2 and the excitonic Bohr radius for the length 
scale a x — h 2 e/fie 2 , where /i = m oe m h/(m oe + m h)- 
The reduced mass fi is determined by optical masses 
where 3/m OK = 2/m^ + l/m|. 

The density functional theory |l4[ expresses the ground 
state energy in terms of electron and hole densities only 
and when applied to the Hamiltonian ([!]) it gives: 



£[n A (r)] =T[n x (r)]+E xc [n x (r)} 
Qxqx' I dvdv 
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n\(r) is the density of the component A [see Eq. (Q)]. The 
Coulomb energy is split into a direct Hartree term and 



an exchange-correlation term. The exchange-correlation 
energy is considered in the local density approximation 
which can be expressed as follows: 



E xc [nx(r)] w / dre IC [ji A (r)]^n A (r) 



(4) 



e xc is obtained from the ground state energy calculations 
of a uniform neutral plasma with equal electron and hole 
densities. In our case charge neutrality can be locally bro- 
ken due to different electron and hole masses. However, 
because of a strong Coulomb interaction the two densi- 
ties differ only slightly and, therefore, it is reasonable to 
calculate e xc substituting an averaged pair density 



n e -h(r) = i ^ [ n e,",<r( T ) + n h,wA v ) 
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Vashishta and Kalia |0| have shown that the exchange- 
correlation energy of a homogeneous electron-hole liq- 
uid is nearly independent of several band characteristics 
of semiconductors, such as valley degeneracy, effect of 
anisotropy, and electron-hole mass ratio. The reason is 
that the anisotropic contribution from the exchange en- 
ergy cancels out the contribution of the correlation en- 
ergy for these band characteristics. Vashishta and Kalia 
further fitted e xc from their self-consistent calculations to 
a simple analytic expression 
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where a = -4.8316, b = -5.0879, c = 0.0152, d = 3.0426, 
and r s — (3/47ra^n e _^) a is the dimensionless distance 
between carriers. 

The densities in Eq. (||) are related to the solution of 
the Kohn-Sham equations (T^| which are self-consistent 
one-particle Schrodinger equations: 



-fiKx + Vx (nx , r) ipM (r) = exiipxi (r) 
/ EA'9A'«A'(r') , 1 
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The sum of the eigenvalues txi over all occupied states 
(A, i) does not give the energy of the molecule but it can 
be directly related to it by: 

E = \^exi - ^qxqx> / dvdv + 

X,i A, A' J 1 1 

+ dr \e xc [n e - h (r)] - 4> xc [n e -h{v)]}n e - h {r) . (8) 
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An approximate solutions to the Kohn-Sham equations 
can be obtained by using a variational ansatz for and 
minimizing the energy (^) . Alternatively, one can choose 
to solve the Kohn-Sham equations directly imposing a 
self-consistency requirement. We calculate a few simple 
cases in both ways and then resort to a variational ap- 
proach in more complicated situations. 



III. BINDING ENERGY CALCULATION FOR 
MULTIEXCITON MOLECULES 

A. Isotropic and equal electron and hole masses 

For isotropic masses (m| = mj: = m K ), we use the 
spherical symmetry of the ground state wave function 
and the effective potential V(r): 



lp\,nlrn(r) = R\,nl{r)Yi m (6 , (j>) 



(9) 



Defining a e . h = m e /m h and R\ : ni{r) = X\,nl(r)/r, we 
obtain the / = Kohn-Sham equations for electrons and 
holes: 
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Xe,no(r) = e e ,nOXe,no(r) , 

Xh,no(r) = e hin oXh,no(r) ■ (10) 



The two equations differ only in the kinetic energy term 
because of unequal electron and hole masses. In this 
subsection we consider the case of equal masses u e -h = 1. 
Then, we have only one Kohn-Sham equation to solve. 
This equation gives the same density profiles for elec- 
trons and holes and, therefore, the direct Hartree term 
vanishes identically The self-consistent solution of the 
Kohn-Sham equation for the lowest eigenvalue is shown 
in Fig. |l| for a molecule formed of 6 electron-hole pairs. 

For the variational solution of the Kohn-Sham equa- 
tions we take the trial functions to be of 'Fermi-Dirac' 
type: 



Is states: R\ s {r) 



2s states: i?2s(f) 
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Here, r%, pi, r 2 , P2 are variational parameters and ni, n 2 
are normalization constants. The choice of the Is wave- 
function is specially suited for large exciton molecules, 
which are described by a 'droplet' model: the density is 
constant up to r% and then drops to zero in a surface 
layer of width pi. The results of the two methods for 
different molecule sizes are summarized in table |. 

The variational and exact solutions of the Kohn-Sham 
equation are in a very good agreement. The difference 



in energies does not exceed 0.5% for all molecule sizes. 
This agreement is not only achieved in the binding en- 
ergies but also in the wave-functions. The comparison 
of the electron-hole pair densities is plotted in Fig. || for 
a molecule formed of six electron- hole pairs. For small 
r the variational trial functions have a linear behavior 
which explains the small difference at the center of the 
molecule between the exact solution and the variational 
solution. The third column of Table | is a measure for 
the stability of the molecule against dissociation into the 
next smaller molecule and one free exciton. We see that 
only the 5- and the 6-pairs molecules are stable against 
this dissociation. The 6 electron-hole pairs molecule is in 
a state where all single particle states of the lowest shell 
are filled up. This represents the most stable configu- 
ration. The energy gain of 24% of E x per one e-h pair 
in a 6-exciton molecule is much higher than a binding 
energy of an ordinary bi-exciton molecule for nondegen- 
erate bands, which is only 1.7% of E x |l7|Jis|| . If we go to 
larger molecules shell effects appear. The extra energy 
cost which is needed to put a further electron into the 
2s shell favors clearly the 6 electron-hole molecule. Shell 
effects appear also if the radii of different molecules are 
compared. Surprisingly, the radius gets smaller when we 
are filling up the Is shell. This feature was also found by 
Wiinsche and co-workers jl2| . The radius of the molecule 
shows a sharp increase if a further electron-hole pair in 
the 2s shell is added. 

The calculated binding energies have to be compared 
to the ground state energy of the metallic electron-hole 
liquid. The ground state energy i? g s .(r s ) was calcu- 
lated in the RPA approximation in Ref. ^ with i?™ n = 
1.55E X at r s = 0.92. The ansatz (|) gives a close value 
-l.6E x , which is reached at the density r s w 1.0. 
Both these values are lower than the maximum gain from 
formation of a 6-exciton molecule which indicates that in 
the chosen approximation there is no intermediate phase 
of a dilute Bose gas of exciton molecules. 

Although it is difficult to give an exact criterion for the 
applicability of the density functional approach to mul- 
tiexciton molecules, it is expected that the theory works 
best at high densities, i.e. for large molecules. Further, 
these densities have to be compared to the density for 
which the ground state energy per pair for an electron- 
hole liquid is minimal r s — 0.92. We expect that in 
the limit of large N, the density at the center of the 
molecule approaches this value. For a finite number of 
electron- hole pairs the density at the center will be higher 
than r s « 0.9, since we must include a surface tension 
which increases the density. To see if such a picture of 
an electron-hole droplet is correct, we split the binding 
energy Eb into a bulk and a surface term: 



pmin 



Eb — E, 



bulk 



(12) 



Esurf • 

We approximate the density by a uniform spherical den- 
sity of radius rjvf- Then the bulk energy can be deter- 
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mined from the Fig. 2 of Ref. ||. The surface energy is 
proportional to r 
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Table || presents the surface tension S evaluated for dif- 
ferent molecules. If an electron-hole droplet considera- 
tion of an exciton molecule is correct, then S should be 
a constant. From table || we can see that this is true for 
molecules formed of 5-7 electron-hole pairs but it fails 
for smaller molecules. 



B. Effect of electron-to-hole mass ratio 

In this subsection we lift the approximation of equal 
electron and hole masses but continue to replace the ef- 
fective mass tensors by isotropic optical masses. The 
numerical procedure of solution of the coupled K-S equa- 
tions ( |l0[ ) was found to be unstable. In view of the re- 
sults of the previous section, we expect, however, to ob- 
tain accurate results using again 'Fermi-Dirac'-type vari- 
ational wave- functions, but this time with different ad- 
justable parameters for electron and holes. Our calcula- 
tions are limited to the most stable molecule formed of 
six electron-hole pairs. We have varied <j e -h from 1 to 0.1. 
The results are plotted in Fig. || and further details are 
presented shown in table III. The total energy of the 6- 
exciton molecule is lowered only by an amount of 0.23 E x 
or 4% if a e -h varies from 1 to 0.1. For CaBg the ratio 
of electron and hole optical masses is (J e -h = 0.89. For 
this value of a e -h the electron and hole densities will be 
only slightly different and there is nearly no change in 
the binding energy Eb- 

We can check validity of the obtained results by com- 
paring them to the exact analytic criteria which have 
been derived in Ref. jlTj for the problem of a bi- 
exciton molecule. Frequently these criteria are not sat- 
isfied by variational solutions. To check them we eval- 
uate the energy functional (|^) in the ground state and 
take it to be a function of <r e -^: E{a e -h). From the 
form of the Hamiltonian (Q) it can be shown [^7) that 
[dE(a e -h)/da e -h]a e - h =i — and that E(a e -h) has to be 
a concave and monotonically increasing function of a e -h- 
We see in Fig. ^ that these two criteria are satisfied by our 
solution, which further supports validity of the obtained 
results. 



C. Effect of mass anisotropy 

According to the results of the previous subsection the 
binding energy changes only very little for unequal elec- 
tron and hole masses. Therefore, for simplicity we as- 
sume that electron and hole dispersion are described by 



the same anisotropic mass tensors. Otherwise, one has 
to deal with too many variational parameters. As in 
the previous section we limit our calculations to the 6- 
exciton molecule. We define the reduced masses 1 / m = 
l/m^r + 1/m^ and 1/m'l = 1 /ml + 1/toJ[ and the ra- 
tio between the components of the reduced mass tensor 
a a = m^/m", with m ± < m)< or a a < 1. 

The variational trial functions are taken to be again 
of 'Fermi-Dirac'-type, but now we use functions with 
a cylindrical symmetry. In different valleys the wave- 
functions are different. For the first valley they are de- 
fined as 



i ) a 



1 + exp 



y 2 +z 2 
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no is the normalization constant and d±, d|| and p are 
the three variational parameters. The wave-functions for 
particles in the other valleys are obtained by cyclical per- 
mutations of x 7 y and z in analogy with the kinetic energy 
operator K\. This implies that the total density distribu- 
tion is not spherically symmetric but is a superposition of 
three ellipsoidal distributions along the coordinate axes. 

The results of the minimization of the energy func- 
tional (j^) are presented in table ^ for 0.2 < a a < 1. 
Varying a a in this range produces an increase in the 
binding energy of 0.51 E x or 6%. For CaBg the mass 
ratio is 0.255, which corresponds to an energy of 7.81 E x 
and thus to an increase in binding by 5%. Thus, mass 
anisotropy has a somewhat stronger effect on binding en- 
ergies of multiexciton molecules in the hexaborides than 
electron-to-hole mass ratio, though both effects give only 
small corrections to the simplest model with equal and 
isotropic masses. 



IV. EXCITON COMPLEXES BOUND TO 
IMPURITIES 

It is well known that excitons in semiconductors are 
attracted to charged donor or acceptor impurities P-p2| . 
We investigate here the case of a monovalent donor im- 
purity, which models the effect of La 3+ -substitution in 
divalent-metal hexaborides CaB@ and SrBg. We add an 
external Coulomb potential to the Hamiltonian ([!]): 



V = I dr^{(r)V x ^ x (r) 

A 



V, = ^ (15) 
er 



and assume equal and isotropic electron and hole masses. 
(In this case there is an obvious symmetry between 
donors and acceptors.) At low impurity concentrations 
the electrons and holes are completely localized at a 
single impurity. The impurity adds a heavy center to 
the molecule. Therefore we choose the variational wave- 
functions in the hydrogenic form: 
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Re(r) = — ^= exp ( -— ) , 

R h (r) = —L^^ex 1? (--), (16) 

V 3lTr h Th V r hJ 

with electrons and holes occupying s-wave and p-wave 
orbitals, respectively, near a positively charged donor. 
The results for the binding energies of multicxciton com- 
plexes are presented in table V. The second column shows 
the total energy of a complex —E B , which includes one 
extra electron in addition to N electron- hole pairs. All 
complexes with N > 1 are stable against dissociation 
into a next smaller complex and a free exciton, since 
\E B (N + 1) - E B (N)\ > E x . To compare the energy 
gain from a formation of a multiexciton complex to the 
energy gain in a dense electron-hole plasma we need to 
subtract from E B the energy of a donor with a single 
electron, which is approximately 2E X . The third column 
shows that excitons gain the most of energy in a 5-exciton 
complex, which is again consistent with a filled shell ar- 
gument. The energy gain per one exciton in a 5-exciton 
complex exceeds the energy gain in a 6-exciton molecule. 
Hence, upon decreasing the band gap charged impurities 
will work as nucleation centers for electron-hole droplets 
and localized excitons will appear before condensation of 
multiexciton molecules in the bulk of a semiconductor. 

The maximum possible energy gain per exciton in the 
5-exciton complex AE = 1.38E X is still below the energy 
gain in a dense electron-hole plasma l-B™™! = 1.55.E X . 
Effective mass calculations have also been done for diva- 
lent donors (acceptors) with 4 excitons being attracted, 
but we did not find any additional energy lowering for 
them. This result means that the direct first-order tran- 
sition occurs under pressure without an intermediate ex- 
citonic phase. However, the two numbers are now closer 
to each other. In such a case so-called central cell cor- 
rections to the effective mass approximation, which have 
been estimated in Ref. || as A£7 C . C , ~ 0.2E X for a single 
exciton, can be sufficient to increase the energy gain per 
one exciton in a donor complex compared to the energy 
density in a bulk e-h plasma. As a result, an intermedi- 
ate state of localized excitons may appear under pressure. 
These localized exciton complexes bound to charged im- 
purities carry uncompensated spin- 1/2, as suggested by 
the shell scenario. It is an interesting problem to under- 
stand whether such an intermediate phase can be respon- 
sible for the unusual relaxation effects in the hexaborides 
& 

Usually, in lightly doped semiconductors (e.g. P in Si) 
localized donor electrons develop antiferromagnetic cor- 
relations between nearest-neighbor donors. This can be 
understood in the following way. For two donor atoms, 
which appear to be neighbors, the electronic hydrogen- 
type orbitals hybridize forming a lower bonding orbital 
and an upper antibonding orbital, similar to a hydrogen 
molecule. The two extra electrons will go to the non- 



degenerate bonding orbital and the Pauli principle re- 
quires them to have opposite spins. We suggest that the 
opposite sign ferromagnetic correlations will develop in 
a localized excitonic phase of doped hexaborides. Their 
origin is in additional multivalley degeneracy of electrons 
and holes in the hexaborides. The bonding orbitals for 
two nearest-neighbor donors have six-fold total degener- 
acy. We expect that in such a case only 4 excitons will 
be attracted to such a double center complex in order 
to fill the bonding electron orbitals. The total spin of a 
double-impurity complex (0 or 1) comes from a partially 
filled hole shell (two 'holes' in the hole shell) . In the case 
of degenerate orbitals the Hund's rule plays the major 
role and produces parallel alignment of spins of the two 
'holes.' Similar scenario applies also to divalent donors 
or acceptors. (In the case of the hexaborides, Ca- or Sr- 
vacancies can play a role of divalent acceptors.) Localized 
moments produced in this way could be a source of the 
unusually fast NMR-relaxation observed experimentally 
|l3f and lead to significant sample-to-sample variations 
of ferromagnetic moments in Cai-^La^Bg 0. 

Estimation of the Hund's splitting in multiexciton 
complexes is an interesting open question. Note, that 
ordinary semiconductors, like Si or Ge, also have a mul- 
tivalley structure of the conduction band. The degener- 
acy of donor orbitals in these cases is partially lifted by 
anisotropic central cell corrections, which select a nondc- 
generate lowest level and, hence, suppress ferromagnetic 
correlations. This effect might be present in multiexciton 
complexes as well. However, the difference between the 
two cases is in the shells which produce the total spin. 
In the multiexciton complexes bound to a donor, uncom- 
pensated spin is formed in the p-wave hole shell. Wave- 
functions for these states vanish at the origin Eq. ( |l6| ) 
and have a small probability in the central cell. There- 
fore, the states in the hole shell are better described in 
the effective mass theory than states in the electron shell 
and their splitting must be less significant. 

V. CONCLUSIONS 

We have shown that the multivalley degeneracy which 
is present in both the conduction and the valence bands 
of the hexaboride materials leads to a number of pecu- 
liar effects: (i) 6-exciton molecules are stabilized due to 
the shell effect. Their binding energies are comparable 
to the exciton binding energy in contrast to a weak bind- 
ing energy of a bi-exciton molecule. The energetics of 
the multiexciton complexes is rather insensitive to the 
details of the band structure: electron-to-hole mass ratio 
and anisotropy in the effective mass tensors, (ii) Mul- 
tiexciton complexes are attracted to charged impurities. 
They have a larger energy gain per one electron-hole pair 
than a 6-exciton molecule. Therefore, upon decreasing a 
semiconducting gap excitons will first condense at donors 
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and acceptors producing an intermediate phase of local- 
ized excitons. Though, the energy of this phase lies some- 
what higher than the energy gain in a dense electron-hole 
plasma, localized excitons can be further stabilized by 
central cell effects or appear as a metastable phase at 
the first-order transition between a semiconducting and 
a semimetallic state, (iii) Uncompensated spins on mul- 
tiexciton complexes show a tendency towards ferromag- 
netic correlations, which is again a multivallcy effect. 

Our theoretical results indicate that a semiconducting 
phase of divalent hexaborides must have a number of 
interesting physical properties. Further experiments at 
ambient and applied pressure can shed more light on their 
relevance to the physics of hexaborides. 

Financial support for this work was provided by Swiss 
National Fund. 



TABLE I. The binding energies Eb of molecules consist- 
ing of N excitons for isotropic and equal e and h masses. 
The equilibrium radius for a given molecule was estimated by 
tm = J dvrn{v)/ J drn(r). All values are given in excitonic 
units. 
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TABLE II. The surface tension for exciton molecules of 
different sizes. 

_N 7 6 5 4 3 

4ttS 0.75 0.79 0.77 0.66 0.49 
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TABLE III. The binding energy of a 6-exciton molecule 
for different values of the electron-hole mass ratio a e -h- 
r e , th, A E , Xh are the optimal values of the variational pa- 



rameters. 


<Je-h 


-E B 


-E B /6 


r e 


A e 




Xh 


1.0 


7.42 


1.24 


1.40 


0.72 


1.40 


0.72 


0.9 


7.43 


1.24 


1.30 


0.75 


1.34 


0.73 


0.8 


7.43 


1.24 


1.32 


0.75 


1.36 


0.72 


0.7 


7.43 


1.24 


1.28 


0.79 


1.34 


0.71 


0.6 


7.45 


1.24 


1.30 


0.77 


1.30 


0.70 


0.5 


7.46 


1.24 


1.28 


0.76 


1.36 


0.69 


0.4 


7.48 


1.25 


1.28 


0.76 


1.38 


0.67 


0.3 


7.52 


1.25 


1.28 


0.76 


1.38 


0.65 


0.2 


7.57 


1.26 


1.30 


0.74 


1.38 


0.62 


0.1 


7.65 


1.28 


1.28 


0.74 


1.42 


0.58 
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TABLE IV. The binding energy of a 6-exciton molecule for 



different values of the anisotropy parameter 


a a . du , d± and p 


are the values of the variational parameters. 








-E B 


-Eb/6 




d± 
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1.0 


7.42 


1.24 


0.52 


0.52 


1.94 
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7.43 


1.24 
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0.55 


1.78 
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7.44 


1.24 
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7.45 


1.24 
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1.76 


0.6 


7.49 


1.25 


0.43 
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1.76 
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7.54 


1.26 


0.40 
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1.73 
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7.61 


1.27 


0.35 
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1.73 
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7.73 


1.29 


0.30 


0.60 


1.69 


0.255 


7.81 


1.30 


0.27 


0.60 


1.66 


0.2 


7.93 


1.32 


0.23 


0.61 


1.62 


TABLE V. The binding energies for different multiexciton 


complexes 


. r s , Th are the values of the variational parameters 


for which the ener; 


ry is minimized. 
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-(E B + 2)/N 
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FIG. 2. The electron-hole pair densities profiles for numer- 
ical (solid line) and variational (dashed line) solutions of the 
Kohn-Sham equations. The results are for the most stable 
6-exciton molecule. 
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FIG. 3. The binding energy of a 6-exciton molecule versus 
the electron-to-hole mass ratio a e - h . 




FIG. 1. The self-consistent numerical Is solution of the 
Kohn-Sham equation for a molecule formed of six elec- 
tron-hole pairs. The resulting effective potential is presented 
on the bottom panel and the lowest eigenvalue is shown by a 
horizontal line. 
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